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ABSTRACT 

We determine the evolution of the galaxy luminosity function (LF) as a function of 
' spectral type using the Autofib redshift survey, a compendium of over 1700 galaxy 

C*") . redshifts in various magnitude-limited samples spanning 6j=ll. 5-24.0. To carry out this 

analysis we have developed a cross-correlation technique which classifies faint galaxy 
spectra into one of six types based on local galaxy templates. Tests and simulations 
show that this technique yields classifications correct to within one type for more than 
90% of the galaxies in our sample. We have also developed extensions of the step-wise 
maximum likelihood method and the STY parametric method for estimating LFs which 
are applicable to recovering an evolving LF from multiple samples. We compare these 
Qlh' methods to the standard 1/V mstx method used in Paper I. 

Applying these spectral classifications and LF estimators to the Autofib sample, 
^ ■ we find that: (i) The spectra and LF of E/SO galaxies show no appreciable evolution 

| out to at least z^0.5. (ii) Early-type spirals show modest evolution, characterised by 

■ a gradual steepening of the faint end of their LF with redshift. (iii) Out to z^0.5, the 

overall evolution of the galaxy population is dominated by changes seen in late-type 
spirals. The characteristic luminosity (L*) of these galaxies appears to brighten with 
redshift and there are signs of strong density evolution (a rapid increase in <£>*). These 
' effects appear to be luminosity dependent so that the LF steepens at higher redshift. 

These trends are accompanied by a steep increase in the median [Oil] equivalent width, 
implying a rapid increase in the star-formation rate with redshift at fixed luminosity — a 
given star-formation rate is found at higher redshift in galaxies of higher luminosity. We 
find that these conclusions are robust with respect to spectral classification errors and 
the luminosity function estimator. Finally, we briefly discuss the correlations between 
our ground-based dataset and a subset of 60 galaxies for which WFPC-2 images are 
now available from Hubble Space Telescope. 

Key words: cosmology: observations - galaxies: galaxies - evolution, large scale struc- 
ture: galaxies - spectroscopy 
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1 INTRODUCTION 



* Please direct correspondence to: Jeremy Heyl, Lick Observa- 
tory, University of California, Santa Cruz, CA 95064, USA 



The Autofib Red shift Survey was introduced in the first pa- 
per in this series ( Ellis et al. 1996) , Paper I) . It is an extensive 
compilation of over 1700 galaxy redshifts drawn from 53 
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separate magnitude-limited slices spanning 11.5<fej<24.0 
About half of the data has been published by o ur ob- 
serving team d u ring 1985-1994 in various papers (Peter- 
son et hi 1985L iBroadhurst, Ellis & Shanks 1988L ICoUess 



et al. l|990|, 1993, |Glazcbrook et al. 1995a| ), and a further 



1026 new redshifts have been secured in the intermediate 
magnitude range 17<6j<22 using the Autofib fibre posi- 
tioner on the 3.9m Anglo- Australian Telescope. The result- 
ing dataset has the advantage of covering a wide luminos- 
ity range up to redshifts z«0.75, and is thus ideally suited 
to studying the evolution of the galaxy luminosity func- 
tion (LF). Details of the photometric selection, observing 
techniques, spectroscopic analysis and a comprehensive cat- 
alogue will be given in a subsequent paper (Paper III). We 
assume H{)=lQQh km s" 1 Mpc" 1 and go=0.5 throughout. 

The main goal of this paper is the determination of the 
evolution of the luminosity functions (LFs) by galaxy type. 
A number of studies at low redshift have investigated the 
LFs for different morphological types, for different colour 
classes, or for different degrees of star-formation activity (in- 
dicated by [Oil] equivalent widt h) . Among the most exten- 



sive of these studies are those of Binggcli, Sandagc fc Tam- 



mann 
and 



1988 



Lin et al. (1996) 



Loveday et al. (1992), Marzke et al. (1994a) 



This work has been directed toward 
understanding differences in the shape of the LF between 
the morphological types, and differences in the LFs of field 
and cluster galaxies. 

In the last few years the size of deep redshift surveys 
has increased to the point where it has begun to be feasible 
to derive the evolution of the LF with redshift. Initial anal- 



yses were carri ed out by Ealcs^l993) 
shi (19^3j| and [Trcyer fc Silk (1994 " 



Lonsdale & Chok- 



using the surveys of 



Peterson et al. (1985)| , |Broadhurst7Ellis fc Shanks (1988) 



Colless et al. (1990)| and powie et al. (1991)} These studies 
clearly indicated evolution of the LF with redshift (as had al- 
ready been apparent from comparison of the galaxy number 
counts and redshift distributions with non-evolving models) . 
However the size of these surveys was not large enough for 
the authors of the LF analyses to reach a clear concensus 
on the form of this evolution. Two new deep redshift sur- 
veys, the Canada-France Redshift Survey of 730 I-selected 
galaxies (Lilly et al. 1995 ) and our Autofib survey of over 
1700 B-selected galaxies (Ellis et al. 1996; Paper I), have 



increased the sample sizes to the point that it is now possi- 
ble not only to recover the evolution of the overall LF with 
redshift, but also to examine the evolution of the LF as a 
function of galaxy type . 

Lilly et al. (1995) measured redshifts for a sample of 
galaxies with I<22 that reaches out to z~l, obtaining LFs 
as function of both redshift and colour. They find that the 
LF evolution is strongly differential with colour, in the sense 
that the LF of red galaxies evolves little if at all out to z~l 
while the blue population shows substantial evolution. Lilly 
et al.'s work probes to the highest redshifts yet explored, 
however their results are limited by the fact that they are 
working with a single magnitude-limited sample, so that 



there is little overlap in the luminosity range spanned in 
different redshift ranges. Thus their results show different 
parts of the LF at different redshifts — the faint end at low 
redshifts, the knee at the median redshift of the survey, and 
the bright-end cutoff at z~l — making it difficult to form a 
coherent picture of changes in the overall shape of the LF. 

The Autofib survey addresses this problem, inherent in 
magnitude-limited samples, albeit at the expense of a lower 
mean redshift than that of CFRS, by combining several sur- 
veys with different magnitude limits, so that we cover a 
broad range of luminosities at all redshifts out to z~0.75 
(see Figure 7 of Paper I). Another difference between this 
work and that of Lilly et al. is that whereas their sample is 
selected in the I band, and so emphasises galaxies dominated 
by older stellar populations, ours is B-selected, and so em- 
phasises the galaxies with ongoing star-formation. Finally, 
whereas Lilly et al. use colour as the indicator of galaxy 
type (dividing their sample into red and blue classes), and 
whereas in Paper I we used [Oil] equivalent width (com- 
paring the LFs for galaxies with EWs greater and less than 
20A), in this paper we use well-defined spectral types deter- 
mined directly from our survey observations. 

The results from the Autofib survey, whilst confirming 
the basic evolutionary trends identified by Lilly et al, are 
more suggestive of the possibility of two distinct popula- 
tions. In particular, the evolving shape of the LF is more 
clearly seen, although to somewhat lower redshift than for 
CFRS, and the trend appears to be largely due to an increase 
in the volume density of 0.1 - 1 L* galaxies with strong [O 
II] emission. In Paper I, the LF was primarily explored as a 
function of redshift for the entire galaxy population. Here, 
we will utilise the spectral classifications to describe the evo- 
lution of galaxy luminosity function in greater detail. This is 
valuable because the spectra themselves offer a glimpse into 
details of the evolving components that cannot be gleaned 
from photometric and redshift data alone. As well as using 
the spectra to isolate different subsets of the data, we will 
coadd the spectra according to various selection criteria in 
order to see if the short-term star formation histories of dis- 
tant galaxies, as revealed by diagnostic absorption lines, are 
different from those of their local contemporaries. This is 
appropriate because the Autofib spectra were obtained at 
unusually good spectral resolution for a faint sample (4-8 
A for bj <22 and 12 A for 22< bj <24 c.f. ~25A for the 
CFRS survey). For a limited subset of our fainter galaxies, 
we also briefly discuss Hubble Space Telescope (HST) im- 
ages and correlations between morphology and photometric 
classifications. 

Our secondary goal in this paper is to discuss these 
spectral classification techniques and also the extensions to 
LF estimation methods we have developed in order to deal 
with the very general nature of the survey sample and the 
extra degree of freedom represented by LFs which evolve 
with redshift. In Paper I these techniques were only briefly 
mentioned in discussing new results on both the faint end 
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of the local LF and the evolution of the overall galaxy LF 
since redshifts z~0.75. 

A plan of the paper follows. The spectral classification 
technique and the tests and simulations carried out to verify 
its precision and reliability are described in §|2| Here we also 
discuss the correlation between photometric properties and 
HST morphologies. In §^ we review existing methods of LF 
estimation and derive extensions to two methods which al- 
low both a clustering-insensitive direct estimate and a best- 
fit parametric model to be recovered for an evolving LF. 
These methods are tested and compared to the standard 
1/Knax LF estimator in In §[| we use the spectral classi- 
fication and LF estimation techniques to determine the LFs 
as a function of both redshift and spectral type. These re- 
sults, and their implications for the physical processes driv- 
ing galaxy evolution, are discussed in Our conclusions 
are presented in §[?]. 



2 SPECTRAL CLASSIFICATION 

In Paper I we stressed the importance of reliable k- 
corrections for the determination of the luminosity function 
over a wide redshift range. The lack of morphological clas- 
sifications or colours for most of the surveys combined here 
means that traditional methods for estimating fc-corrections 
cannot be uniformly applied. Only the DARS survey has 



been classified morphologically and only Lhe LDSS-1 and 



LDSS-2 samples have b 3 — i j colours. Applying a mean k- 

correction or defining a new passband probably would not 
be fruitful as the galaxies span z~0-l, so an extremely blue 
or red galaxy at high redshift would have its luminosity in- 
correctly estimated by over a magnitude, and using a mean- 
redshift could lead to errors nearly as large. Furthermore, 
the volume weighting implied would be even more uncer- 
tain. As a major goal of our survey is to estimate the LF 
as a function of spectral class, we need an approach that is 
both uniformly applicable and more reliable. 

In principle, a fc-correction could be 'read' directly from 
each flux-calibrated spectrum. In practice, the calibrations 



are toolUriceftain to allow the continuum shape of the spec- 



tra to bre matched Lu a standard GED, as is demuusLiaLed by 

a comparison of synthetic colours estimated from the spec- 
tra with those obtained from direct imaging. However, on 
smaller scales where spectral features can be measured, the 
data are considerably more reliable. Accordingly, we have 
developed a technique for classifying each galaxy spectrum 
in the survey which succeeds very well for but the most noisy 
spectra in the sample (those from the LDSS-1 and LDSS-2 
surveys, for which we used the bj — rf co lour to estimate 
the fc-correction as in Colless et al. (1990)). This technique 



was briefly outlined in Paper I (§3.1); here it is described in 
detail and validated by a variety of tests. 



2.1 The Cross-Correlation Method 

To determine the spectral classification of each of the galax- 
ies, we chose to cross-correlate the survey spectra against 
those of the Kennicutt (1992a, 1992b) spectral library. This 
library is highly appropriate for this task, not only in spec- 
tral resolution and wavelength range but also because the 
spectra represent the light integrated over a large portion 
of the local galaxies rather than of the central regions as is 
usually the case. The physical apertures involved are closely 
matched to those of Autofib and LDSS at their mean red- 
shifts. The 2-arcsec effective aperture of these surveys corre- 
sponds to a physical size of 2-8 h' 1 kpc for redshifts between 
0.1 and 1.1, well matched to the physical apertures used by 
Kennicutt (1992a, 1992b). 

The cross-correlation technique works as follows. First 
both the Kennicutt template and survey spectra are 
smoothed on a scale of 100A in the observer's frame and 
these smoothed versions are subtracted yielding continuum- 
subtracted spectra. The spectra are then rebinned to 2A per 
bin in the rest frame of the galaxy. The cross-correlation is 
defined to be 



(Ei ti°i 



(1) 



where ti and Oi are flux values in the 2A bins in the template 
and observed spectra respectively. This procedure is similar 
to the method der ived independently by Zaritsky, Zabludofi 



Willick (1995). The survey spectrum is classed accord- 



ing to the morphological type of the template with which 
it most strongly correlates. Each Kennicutt template has 
an associated morphologi cal classification give n in Kenni- 



cutt (1992b). We find the [King fc Ellis (1985)| fc-correction 
that most closely corresponds to the Kennicutt classifica- 
tion, thus determining the fc-correction for each galaxy over 
a range of redshifts for which it might be detected within 
the particular sub-survey given the apparent magnitude lim- 
its. We cannot use the Kennicutt data to directly calculate 
the bj fc-corrections for z > 0.05, as they do not extend into 
the ultraviolet below 3650A. On the other hand, the spectral 



energy dis tributions (SEDs) of Pence (1976) and King & El- 



lis (1985), which do extend into the UV, lack the small-scale 
spectral features that make t he cross-correlation method 
possible. Kinney et al. (1993) have published a UV spec- 
troscopic atlas of star-forming galaxies; unfortunately, the 
uncertainties in these data are too large in the wavelength 
range critical to estimating fc-corrections (2800-3700A) for 
this atlas to be useful for our purposes. Additionally, the 
IUE aperature is not well-matched to the distant popula- 
tion. By using the Kennicutt spectra as templates and the 
appropriate SEDs for the matched types, the fc-correction 
can be readily found. 

If one assumes that the noise level of the survey spec- 
tra is uniform with wavelength (a good assumption for the 
faintest spectra), the cross-correlation coefficients may be 
directly converted to a value of \ 2 '■ 
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Table 1. Success rates of the classification algorithm. 



Fluxed Unfluxed 



S/N 


Template 


Type 


Template 


Type 




Correct 


Correct 


Correct 


Correct 


4.05 


100% 


100% 


81% 


84% 


2.03 


94% 


94% 


76% 


83% 


1.22 


88% 


90% 


62% 


72% 


0.82 


82% 


86% 


52% 


68% 



X 2 = ^E°^ r ) ( 2 ) 

i 

Clearly the best template match will yield a minimum \ 2 
value. In principle the x 2 estimator enables the method to 
be made more sophisticated by including both a distribution 
of prior probabilities among the template galaxies and a list 
of those templates which are a near-equal match. This in- 
formation may provide the basis of an error estimate for the 
fc-correction of each galaxy, but is beyond the scope of this 
analysis, for which we used the 'raw' r-values to determine 
the best-matching galaxy. 




Figure 1. Input versus output classifications. The height of the 
surface represents the fraction of objects within each input class 
that were classified into each output. Red ellipticals and early 
spirals can be confused with blue ellipticals; and intermediate 
spirals can be confused both with earlier and later types. 



2.2 Simulated Tests 

To verify the algorithm, we performed two series of tests. 
The first involved simulating the procedure using fluxed 
spectra. We randomly selected a Kennicutt spectrum, nor- 
malised it to a particular mean number of counts per bin 
and added a sky signal. Using the total count per bin, a ran- 
dom Gaussian deviate was chosen and this noise added prior 
to sky-subtraction. By repeating this process 100 times, we 
created an ensemble of test spectra of known signal/noise 
ratio in each pixel. Each of these spectra was processed as 
above and the success rate in returning the correct Kenni- 
cutt class was calculated. As Table |l| shows, the success rate 
is encouraging. 

Next we tried to simulate the unfluxed spectra in the 
catalogue and also took into account the effect of varying the 
redshifts of the simulated spectra. Test spectra were gener- 
ated as before except that each spectrum was multiplied 
by a response function which is zero outside 3600-7200A 
and increases quadratically by a factor of two from the edge 
to the central wavelength. This response function can be 
blueshifted by a factor of x 1-1.6 equivalent to a redshift of 
0-0.6. Table 1 shows the success rates for these tests. 

Figure h examines the results of the "unfluxed" simu- 
lations in more detail. The most striking feature of the dis- 
tribution is the diagonal ridge line which traces the correct 
classifications. Classes #2 (red ellipticals) and #4 (early spi- 
rals) are sometimes confused with class #3 (blue ellipticals). 
Furthermore, class #5 (intermediate spirals) are sometimes 
classified as class #3 (blue ellipticals) and as class #7 (star- 
bursting spirals). However, the overall agreement is excel- 
lent. 



As mentioned above, each fc-correction class corre- 
sponds to several Kennicutt spectra, each with varying 
strengths of spectral features. Consequently, with the ad- 
dition of noise, the spectra can be confused across a k- 
correction class. In the case of spectra for intermediate spi- 
rals which contain weak features of both earlier and later 
classes, the classification may be significantly in error. 

Finally, we performed the blueshifting described above 
to simulate how observing galaxies at various redshifts 
through a fixed wavelength range can affect the success 
rate. Figure |^ shows that the error rate increases because 
a decreasing portion of a galaxy's spectrum overlaps with 
the Kennicutt templates. We can only expect a satisfactory 
classification to z ~ 0.6. For higher redshift galaxies in the 
LDSS-1 survey, and for all the LDSS-2 survey galaxies, we 
use colour to determine spectral class. 



2.3 Sample Spectra 

Figures |^ and ^ illustrate examples of the spectral classifica- 
tion technique. In each case the lower curve is the observed 
spectrum and the middle curve is the Kennicutt spectrum 
selected as the closest match. The continuum has been sub- 
tracted from both spectra as described earlier and the ob- 
served spectrum has been smoothed further to accentuate 
the features. The upper curve is the product of the two spec- 
tra smoothed over 20 bins to show which features contribute 
most strongly to the total cross-correlation. 

Examining the brighter survey galaxies (Figure ^|) there 
is strong correspondence between the observed spectra and 
the best templates. Intermediate-type galaxies are the most 
challenging as seen from the middle panel of Figure ^. How- 
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Figure 2. Classification error rate as a function of redshift. The 
error increases markedly with redshift, as the identifying features 
of the spectra are lost off the red end of our simulated spectro- 
graph. The horizontal line is the mean error rate of the various 
tests: 29%. 



ever even with faint galaxies (Figure ||) the algorithm per- 
forms well. However, the limited overlap at higher redshifts 
is a stumbling block (as predicted by Figure b|). Presently 
there is no equivalent spectroscopic atlas covering near- 
ultraviolet wavelengths, and in any case there are few strong 
features blueward of the [Oil] emission line that would be 
helpful. For more distant sources, the most promising option 
would be to extend the observing window redwards utilis- 
ing night sky suppression techniques to limit the deleterious 
effects of the OH emission. 



2.4 Tests with Real Spectra 

A subset of the galaxies in the catalogue have bj — rp colours 
and this provide an independe nt test of the classi fication al- 
gorithm. The LDSS-1 galaxies (dolless et al. 1990) lie within 
21< B <22.5 and are amongst the faintest in the survey. 

We use a two-stage test: the correlation of observed 
colour against fc-correction class and correlation of 'pre- 
dicted' colour against observed colour. We calculate the 
b.j - tf colour in the observer's frame by using the classifi- 
cation of each LDSS spectrum using our method and then 
determining the colour of the appropriate SED at the red- 
shift of the galaxy. Figure |E] shows the results of both these 
tests. 

The trend of colour versus class (first panel) exhibits a 
large spread which could arise from observing errors, colour 
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Figure 3. Classifying bright spectra. Three continuum- 
subtracted spectra selected from the bright section of the Aut- 
ofib survey — from top to bottom, an early-type galaxy (B=19.67, 
2=0.127), an intermediate-type galaxy (B=17.09, 2=0.047) and 
a late-type galaxy (B=19.57, 2=0.049). The lower curve in each 
panel traces the the observed spectrum after the continuum 
has been subtracted, best-fitting (continuum-subtracted) spec- 
trum from the Kennicutt atlas. The middle curve is best-fitting 
(continuum-subtracted) spectrum from the Kennicutt atlas, and 
the upper is the cross-correlation contribution as a function of 
wavelength. 
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Continuum-Subtracted Early-Type Spectrum; B = 2 1 -68. z = 0.45j 




2500 3000 3500 4000 4500 5000 

Wavelength (A) 



Continuum-Subtracted Late-Type Spectrum; B— 21.87, z = 0.564 
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2500 3000 3500 4000 4500 

Wavelength (A) 

Figure 4. Classifying faint spectra. Two spectra selected from 
the faint section of the Autofib survey. The upper panel is an 
early-type galaxy (B=21.87, z=0.564). The lower panel is a late- 
type galaxy (B=21.68, 2=0.453). The arrangement is as in the 
previous figure. 



corrections and misclassifications. The second panel shows 
the predicted colour against the observed colour thereby re- 
moving the effect of redshift. Although the trend is remark- 
ably good the mean-absolute error is approximately 0.4 mag. 
Budgeting the errors for this multistage process is difficult, 
as they arise in several places. The bj and magnitudes 
each have on average an uncertainty of 0.1 mag leading to a 
total error in the colour of 0.15 mag. Furthermore, a few 
galaxies have colours bluer than any of the fc-correction 
SEDs; this is worrisome, but the effect is small and limited 
to only a few galaxies. 

The remaining spread of approximately 0.3 mag may be 
due to either the variance of the i ntrinsic colours of ga lax- 
ies within a given spectral class ( deVaucouleurs 1977) or 



* i * 




1 2 3 

bj — r F (observed) 



Figure 5. Results for LDSS-1 spectra. The upper panel depicts 
the spread in observed bj — rp colour for the galaxies in the 
LDSS-1 survey against the fc-correction classification. Each small 
square represents a single galaxy. The large circles show the trend 
in the median colour for each class. The lower panel shows colour 
predicted from the fc-correction classification at the redshift of 
each galaxy against the observed colour of the galaxy. Each square 
represents a galaxy in the survey. 
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to misclassifications. Although it depends on fc-correction 
class, the error in the fc-correction is approximately 1.5 times 
greater than the colour error, so if all the spread in this plot 
were due to misclassifications, the fc-corrections for these 
galaxies would have an RMS error of 0.6 mag. The LDSS-1 
data are the most noisy of the survey and we expect some 
scatter in the intrinsic galaxy colours; consequently, we ex- 
pect the mean misclassification rate and the fc-correction 
error to be much lower over the entire survey and an unim- 
portant source of error in the luminosity function estimates. 



2.5 Testing with HST 

Given the limited usefulness of comparing spectral classi- 
fications with those based on broad band colours, it is in- 
teresting to compare both classifications with those based 
on morphologies derived from the Hubble Space Telescope 



(HST) images. For the brighte r redshift surveys (Peterson 



et al. 1985, Loveday et al. 1992) fc-corrections have normally 
been assigned on the basis of a Hubble type determined from 
visual inspection and thus it is important, where possible, to 
examine the correlation between Hubble type and the diag- 
nostics we have employed in this work to classify our spectra 
and faint galaxies. 

As part of a major programme to morphologically re- 
solve the faint field population, we have imaged about 80 
galaxies in the fields sampled by our survey with 21 < 
bj <24 with the HST WFPC-2. Typical integration times 
were 4-6 orbits in the F814W filter. Of these 80, 64 have 
spectra and redshifts as part of the Autofib survey. Most 
of these are drawn from the LDSS-2 catalogue although a 
few LDSS-1 and fibre sample galaxies were also observed. 
At these apparent magnitudes the spectra are generally too 
faint or highly redshifted for the cross-correlation technique 
to be a reliable classifier (c.f. Figure 2) and the bulk of the 
classifications are thus based on bj-rF colours. 

Each of the 64 galaxies has been assigned a morpho- 
logical type by one of us [RSE] according to the mor- 
phological sch eme introduced for analy ses of the Medium 
Deep Survey (Slazcbrook et al. 1995b ) and Hubble Deep 
Field (Abraham et al. 1996). Briefly, galaxies were classi- 
fied visually on the F814W image according to a scheme 
whose internal resolution adapts to the varying quality 
and size of the individual images. Although the scheme 
has 12 one-dimensional classes, in this comparison we 
have collated the data into 5 broader classes, namely stel- 
lar/compact, spheroidal, early /intermediate spiral, late- type 
spiral/Irregular, merger/unclassed. In the following we ex- 
amine how well these HST morphologies correlate with the 
classes used in the rest of the paper noting that these are, 
largely, based on colour. 

Figure ^ shows the bj-rp colour-redshift relation with 
each datapoint denoted according to its HST morphology, as 
well as a direct comparison of the HST morphological type 
with the fc-correction class based on the broad-band colours 



or spectroscopic cross-correlation. In these plots 3 QSOs in 
the HST sample have been removed but objects which defied 
classification in one or both schemes have been included. 

As before, recognising that morphology 
and colour/spectral class do not perfectly align even with 
high signal/noise local data, we can consider the agreement 
to be satisfactory; 52% of the HST-classed galaxies lie on the 
ridge line, sharing identical classifications. The largest sys- 
tematic deviation appears to be with the morphologically- 
defined spheroidal galaxies where a significant proportion 
have spectral or photometric properties which class them as 
early-type spirals. On closer inspection these sources include 
galaxies which appear to lie at the boundary of the SO and 
early-type spiral region. However, as discussed later, there 
do appear to be a number of cases where HST-classed ellip- 
ticals have slightly bluer colours than expected and this is 
accompanied by [O II] emission, so the blurring of classes in 
the spheroidal/early- type spiral region is partially a physical 
effect. 

More extreme misclassifications than those which trans- 
fer data to adjacent bins are rarer and, upon inspection, 
easier to understand. For example, there are quite a few 
compact objects which are HH-galaxies with extremely blue 
colours and spectra indicative of late-type spirals. Under our 
simple morphological scheme, these are binned with the el- 
lipticals since this classification boundary for distant objects 
is difficult to discern. There is also one compact red object 
with an emission line. 



3 LUMINOSITY FUNCTION ESTIMATORS 

Once the luminosities of survey galaxies have been deter- 
mined (see § 3.2 of Paper I), the next logical step in our 
analysis is to construct the luminosity function as a func- 
tion of class. A number of possible LF est imators have been 



discussed in the literature (see rev iews by Felten (1976) 



stathiou, Ellis & Peterson (1988)). In the particular case of 
the Autofib survey, we are interested in an estimator which 
is well-suited to the deep pencil-beam nature of the survey 
which may imply greater large scale structure fluctuations 
than would be the case for panoramic all-sky surveys. 

In Paper I, we presented LF results based on a sim- 
ple estimator, the 1/Knax method ( Schmidt 196£ ), arguing 
that although the redshift distribution in individual pencil 
beams may suggest large departures from that expected for 
a homogeneous distribution, the large number of individual 
beams (53 in 9 widely distributed regions) should minimise 
such difficulties. Here we will explore this assumption more 
carefully by developing a modified version of the Step- Wise 
Maximum-Likelihood (SWML) method. That method, and 
its derivatives, has become a popular tool to overcome the ef- 
fects of fluctuations by assuming the LF has the same shape 
across all density fluctuations in the survey. Although useful 
in local surveys, as we are interested in examining potential 
evolutionary changes in the LF this restriction is unaccept- 
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able. In the following sections we derive and compare various 
alternate methods for estimating the LF which do not im- 
pose this constraint. 



3.1 The 1/Vmax Method 

We begin by briefly reviewing the 1/V m ax method used in 
Paper I. The first step in this method is the calculation of 
the total volume within which the object could have been 
observed. In a single magnitude-limited survey (see § 3.2 of 
Paper I) : 



V m 



H 



A 



d L dz 



(1 + z) VI + ' 2 Qoz 



(3) 



where z m i n and z max are the minimum and maximum red- 
shifts from which the object could have been observed in the 
survey considering the distance modulus and fe-correction. A 
is the area of sky surveyed in steradians. In its simplest form, 
the luminosity function is obtained by collecting the sources 
in bins of constant magnitude and summing the (Vmax) 1 
values in each group. 

Avm & hiahcall (lysu/ describe how to combine more 
than one sample coherently in a V/ Vmax analysis. The new 
variable, denoted V a is simply the sum of V max over all the 
surve ys in which t he object could have been observed. Fi- 
nally, Eales (1993) describes how to use this variable to con- 
struct the luminosity function as a function of redshift. The 
analysis proceeds according to Schmidt's method with the 
exception that the data is also binned in redshift. Here z m in 
is the minimum redshift at which the object could have been 
observed in the magnitude-limited sample and be in the red- 
shift range of interest. The method works similarly for the 
maximum redshift. 

As mentioned, these variants of the traditional 1/Vmax 
estimator are each affected by the presence of clustering in 
the sample. Clustering will lead to a poor estimate of the 
faint-end slope of the luminosity function. 



Star/Compact E/SO Early-Mid S Late S Unclass 



HST Classification 



Figure 6. Correlations with HST morphology. The upper panel 
shows the observed bj — rp colour-redshift relation for 55 galax- 
ies in the LDSS-1 and LDSS-2 surveys. Different symbols have 
been used according to the 5 HST morphologies discussed in the 
text (see inset). The curves represent the expected trends for the 
adopted set of fc-corrections. The lower panel shows a direct com- 
parison between HST morphology and photometric/spectroscopic 
class including those sources which defied classification in one or 
both systems. The area of each circle is proportional to the sample 
size in each category and is indicated numerically. 



3.2 Clustering-Insensitive Methods 

Several authors have introduced estimators, which although 
slightly biased, reduce the effects of clustering on the result- 
ing LFs. Each of these techniques attempts to maximise the 
likelihood of observing a particular set of galaxies by varying 
parameters of the luminosity and completeness functions. 
The methods differ in the way the assumed probability of 
observing a particular galaxy is calculated. 

Adopting Poisson statistics is an obvious way of pro- 
ceeding. This assu mption was first app lied to the analysis 
of quasar samples ( Marshall et al. 1983 ) and later rec ast to 
analyse the CfA redshift survey ( Choloniewski 1986 ). The 
probability of observing k galaxies in the interval dMdz in 
an area of sky d£l is 



(4) 



Luminosity Function Evolution by Spectral Type 9 



where 



3.3 Deriving the SSWML Method 



A = ~<f>(M)p(z)dMdz<Kl 
n 



(•») 



and n is the average number density of the survey. By bin- 
ning (j)(M) and p(z), these two functions may be estimated 
without assuming particular forms. The only task that re- 
mains is to maximise the total probability of the survey (the 
product of all the Pk in each of the bins) while varying the 

values of 4>(M) and p(z). 

The C-method (Lynden-Bell 1971) adopts a different 



approach. By considering the plane of redshift versus abso- 
lute magnitude, this method uses the fact that the ratio of 
the number of galaxies observed between L and L+dL to the 
number brighter than L is proportional to the ratio of the 
number of galaxies actually in the field in this range to those 
brighter than L multiplied by a weighting factor to account 
for the differing volumes sampled. The method generates a 
cumulative luminosity function without normalisation. The 
differential luminosity function may be derived by fitting an 
appropriate model. 



The STY method ( |5andage, Tammann fc Yahil 197S| ) 
twists this plane around, ft examines the probability that a 
galaxy observed at a redshift zq is brighter than M: 



P(M, z ) 



<KM')D(zo)f(m')dM' 
<p(M')D(z )f(m')dM' 



(6) 



where f(m) is the completeness of the survey at apparent 
magnitude m and D(z) is the density of galaxies at red- 
shift z divided by the mean density of galaxies (p). Taking 
the derivative of this equation with respect to M yields the 
probability density for finding a galaxy with absolute mag- 
nitude M in a magnitude-limited survey. This probability is 
directly proportional to the density of galaxies with that ap- 
parent magnitude and inversely proportional to that which 
could have been observed at that particular redshift, 



p k oc <p(M k 




M bright( 2 fc) 



cj)(M')dM' 



(7) 



M faint( z fc) 



Here we have replaced the function f(m') with a function 
that is zero outside the magnitude limits of the survey and 
unity within (i.e. assuming that the survey is 100% com- 
plete). As the redshift of the galaxy is fixed, the discontinu- 
ities in f(mf) correspond with the range of absolute magni- 
tudes beyond which no galaxies at this redshift could have 
been observed. 

Finally, the Step- Wise Maximum-Likelihood metho d 
follows from equation | ( fcfstathiou, Ellis fc Peterson 1988|) . 
Our modification of the SWML method (SSWML) comple- 
ments Choloniewski's method, but for a different statistical 
model. 



The aforementioned clustering-insensitive methods by de- 
sign probe the LF as a function of luminosity only. To un- 
derstand the evolution of the LF with redshift, this restric- 
tion must be removed. Here we will derive two new methods 
that avoid this restriction. They are based on, respectively, 
the STY and SWML methods. Our generalisations of the 
locally- valid estimators will be denoted with a prefixed 'S' 
for spatial, i.e. SSTY and SSWML respectively. 

In the derivation, it becomes apparent that these gen- 
eralisations provide some additional rewards. For example, 
they provide a straightforward prescription to combine var- 
ious surveys coherently and to determine the absolute nor- 
malisation. 

The derivation of the SSWML method begins with 
equation ^ of the STY method, modified by two generali- 
sations: 

• p(z, M) 7^ pD(z)(j>(M) (i.e. luminosity and density evo- 
lution are separately permitted), and 

• f( m ) — fi(m), where Q(m) is the area of sky sampled 
at apparent magnitude m, accounting for sampling rate and 
mean completeness. 

These two generalisations will allow the determination o f the 
LF as a function of redshift (as in Choloniewski 1986) and 
the use of many surveys in a single coherent determination 
of the luminosity function. 

The relevant probability is that of finding a galaxy in 
the survey brighter than M and closer than z: 



P(M,z) = 



mlp(^M'Mm')^dM'dz' 
Io°° JZ P(*'> M 'Mm')%dM>dz' 



(8) 



Again, taking the derivative of this equation yields the prob- 
ability density: 



Pk 



d 2 P{M,Z) 
dMdz 

oc p(zk,Mk) / 



(9) 



oo />oo 



J -oo 



r (z',M')Q(m')^dM'dz' 



where fi(m) is the solid angle sampled at the apparent mag- 
nitude m, and m is the apparent magnitude corresponding 
to M and z considering the distance modulus and effects 
of the fc-correction. Consequently m may best be written 
m(z,M;Ck) where Ck is the fc-correction class of galaxy fc. 
This equation forms the basis of the SSTY method. 

Figure [?] compares Equations |^ and [] on the absolute 
magnitude versus redshift plane. The lower double integral 
is simply the number of galaxies that one would expect to 
observe in a combined survey given a LF. If the function 
Q(m) is simply a series of steps (as in Figure (?]), there is 
only one fc-correction class, and the trial LF is integrable at 
a given redshift (like an evolving Schechter function), the 
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° Unobserved Galaxy 



• Observed Galaxy 




stant everywhere, we can choose to perform the latter inte- 
gral at the redshift of the observed galaxy (zk), yielding 



Pk oc n(zk)<j>(Mk) 



(13) 



rlV f M bright( z k) 

n{z)Q,—dz \ <j}{M)dM 



which reduces to the STY result. 



We now proceed in the spirit of the S WML ( Efstathiou 



Ellis & Peterson 1988) method and assume that p{z,M) is 



Redshift *> 

Figure 7. Geometric comparison of the STY and SSTY methods. 
In the STY method, the probability of observing a galaxy is the 
ratio of the number of galaxies at the galaxy's absolute magnitude 
(the dark-shaded region) to the total number that could have been 
observed within the magnitude limits at the same redshift as the 
observed galaxy (the light-shaded region). The extension to the 
SSTY technique is straightforward. The probability of observing 
a galaxy at a particular redshift is the ratio of the density of 
galaxies at the observed redshift and absolute magnitude (again 
the dark-shaded region) to the total number of galaxies that could 
have been observed at any redshift (the region enclosed by the 
magnitude limits). 



integral may be most efficiently calculated as 



xrpred 
JV gal 



> Q,j / ——dm'dz x 

~ Jb. ■ dz 

j 3 ,min 

p(z , m' — dmoduius(z) — k(z)) 



(10) 



where the sum is over the subsurveys and the jih survey 
samples Qj steradians in an apparent magnitude range from 
-Bj.min to Bj |max . Although the function p(z,M) may not be 
separable, if the evolution of the LF is posed as an evolution 
of separable Schechter parameters, p(z, M) takes the form 
of an integrable Schechter function at all redshifts. 

Although now complex, if we reinstate the assumption 
that the density is separable and that we have a single 
magnitude-limited sample, this formula reduces to Equa- 
tion [jj 

/roc r-oo 
J n ^~d7 J <M M M m ) dM ' d2 ' (ii) 

We can calculate fi(m) for z k 

0(m) = / if A-if a i nt (Zfc) < M < Afbright(zfc) ^ 

1 otherwise 
Since we seek an estimator for 4>(M) which is assumed con- 



defined for a two-dimensional array of steps in both redshift 
and luminosity, 



p(z, M) = 2_^ W{z ~ Zi,M- Mj)pij 



(14) 



if 



Az/2 < z < Az/2 and 
-AM/2 < M < AM/2 . (15) 
otherwise 



where 
W(z,M) = 

Substituting this relation into the formula for pk gives 

p k oc ^ W(z k - Zi, M k - Mj)pij / (16) 



+Az/2 /-Mj+AM/2 



-Az/2 J Mj-AM/2 

Next we calculate the logarithm of p k 

lnp fc = >J W(z k - Zi,M k - Mj) In p tj - In I k 

ij 



Q.(m)^dMdz. 

dz 



(17) 



where we have replaced the denominator with I k - The sub- 
script k is retained because the value of the integral de- 
pends on the fc-correction for the particular galaxy observed 
through f2(m(z, M; c k )). 

The likelihood of observing the entire survey is the 
product of the likelihoods for each galaxy 



n 



p k , 



and taking the logarithm, 



InPsurvoy = ^ lnpfc. 



(18) 



(19) 



Because we would like to maximise the survey likelihood, 
we examine the derivative with respect to p pq and search 
for stationary points: 

dlnP sm . vcy _ ^ d\np k 



dppq ^ dppq 

k 



o. 



Differentiating Equation [It], 
dlnpfc _ W{z k - Zp, M k - Mg) 



dp p 



Pp'i 



1 dl k 

I k dp pq 



(20) 



(21) 
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Referring to Equation [l6| we see that 

, r r z p + Az/2 r M q +AM/2 

f±- = / / n(m)^dMdz 

"■PPI Jz p -Az/2 JM q -AM/2 dz 

and 

, _ dlk 

1 k — > Pp9 j 

^ dp pq 



(22) 



(23) 



Summing over the derivatives (Equation ^lj) gives 

dlnP survey \ / 1 . . 1 dlk 

= > W(z k - z p , M k - M q ) - — 



.(24) 



and setting this to zero yields 

which is familiar in form to the SWML result, 
</> q AM = ^ W ( M * ~ M i 

k 

< H[M g - M iaint{zk) ] 

k I , 



(25) 



(26) 



jAMH[Mj - Mi 



faint(z fe )J 



since we can identify with the sum over j in the denom- 
inator and dl /dp P q with H[M q — Mf a i nt ( z .)] also in the de- 
nominator. 

Calculating the values in the denominator of Equa- 
tion [25] is less straightforward than in the SWML method, 
but Equations ^ and ^ show that this denominator is a 
function of the p pq , the fc-correction class of galaxy k, the 
cosmology (through dV/dz and the distance modulus), and 
the details of the survey (through f2(m)). Therefore we can 
calculate the values of the dlk/dp pq for each fc-correction 
class before beginning the iterative solution to Equation |25| 

Furthermore, since the integrals I k are over volume, it 
is straightforward to calculate the total number of galaxies 
expected given the current values of pij, 



prod _ V"^ r dlk 



(27) 



ijk 



where fk is the fraction of the galaxies observed in each 
fc-correction class. All p pq can then be normalised so that 
Equation yields the number of catalogue galaxies. The 
normalisation can be done either at each iteration by multi- 
plying the right-hand side of Equation || by n^ CIVcd /™llf , 
or by multiplying this ratio after the final iteration. The al- 
gorithm converges more quickly (two or three iterations) if 
the normalisation is performed at each step. 

Errors may be estimated following procedures discussed 



by [Baunders et al. 199(j : 

f(log/9 Pq ) = (lnlO) _1 cr(lnppq) 



= (lnior 1 



fr In P 

u 111 1 survc 

{d\npp q y 



-1/2 



(lnlO) -1 x 



Y \ W(z k -z p ,Mk-M q 

I- ^ 



(28) 
(29) 

2 \ I - V2 



dlk I \ ^ dl k 

ppq d Ppq / 2^ pij d Pij 



Additionally, it is straightforward to estimate upper limits 
for bins in which no galaxies are observed. 



1/2 



J2k [ dI k/dp Pq jYLn pijdh/dpij. j 



(30) 



In this derivation, Q(m) includes several important fac- 
tors. The simplest way to define Q(m) is to calculate the area 
surveyed in each subcatalogue and then multiply this area by 
the sampling rate and completeness. f2(m) then represents 
a series of steps with a jump at the bright and faint limits of 
each subcatalogue. However, we can generalise this in two 
ways. Firstly, since we explicitly calculated dlk /dpij for each 
galaxy type, we can introduce completeness as a function 
of fc-correction type by defining various fifc(m). Secondly, 
we can calculate the completeness rate within a subcata- 
logue as a function of apparent magnitude and account for 
this in obtaining fi(m) or flfc(m). These two generalisations 
make this technique an extremely versatile tool in analysing 
galaxy catalogues. 



4 TESTS AND COMPARISONS 

To test the 1/Vmax and SSWML methods in a large sam- 
ple like the composite DARS, AUTOFIB and LDSS-1/2 
surveys, we generat ed a random galaxy catalogue from a 
Schechter function (Loveday et al. 1992). In the first test 
1800 galaxies were selected in a standard cosmology, 300 
from each of the following magnitude ranges: 11.0-17.3, 
17.0-19.7, 19.7-20.5, 20.3-20.8, 20.8-22.5, and 22.5-24.0. All 
the galaxies were assigned zero fc-correction and the density 
of galaxies was doubled beyond a redshift of 0.2 to crudely 
simulate density evolution. 

Figure [] illustrates the results of this simulation. It can 
be seen that the two methods generate nearly identical LFs, 
although the SSWML method produces smoother results. 
The bootstrap errors of the 1/Vmax method agree remark- 
ably well with errors derived for the maximum likelihood 
technique. Although in principle the SSWML technique has 
the additional advantage of determining upper limits which 
further constrain the evolution of the faint-end slope, these 
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Figure 8. Test catalogue of 1800 galaxies. The upper panel de- 
picts the results of the 1/V ma x analysis on the random galaxy 
catalogue and the lower panel shows the results of the SSWML 
method. In each panel the lower curve is the nearby (z < 0.2) 
luminosity function, and the upper is the distant one. The er- 
rorbars for the l/Vmax method are generated using a bootstrap 
technique, while the errors and upper limits for the SSWML are 
determined as described in the text. 



upper limits increase with a slope of a = — 2 and so in prac- 
tice are not very interesting much fainter than the observed 
points. 

To test whether either of the techniques are biased and 
to find an external error measurement, we generated an en- 
semple of ten catalogues derived from the same evolving LF 
as depicted in Figure^. Each catalogue contained 1,800 ob- 
jects sampled as in the previous test. We found that neither 
method is systematically biased relative to the input lumi- 
nosity function and that the internal error estimates satis- 
factorily approximate the scatter of the luminosity functions 
derived for the ensemble. 



4.1 Sensitivity to Clustering 

To test the sensitivity of the algorithms to clustering, three 
additional galaxy catalogues were generated with an over- 
density of galaxies at a redshift of 0.05 with a dispersion in 
redshift of 0.003 in one of the fields in the 17.0-19.7 sur- 
vey. The 'cluster' LFs were assumed identical to that in the 
field, but were normalised such that a set fraction of the 300 
galaxies in the 17.0-19.7 magnitude range were in the clus- 
ter. The three surveys had clustered fractions of 45%, 65% 
and 85%. 

Figure ^| shows the results of this simulation for the most 
populous cluster for a variety of bin widths. The 1/Knax and 
the narrowly-binned SSWML methods both over-predict the 
faint end of the LF. However, as the bin width of the SS- 
WML method is increased, the over-prediction decreases. 
The 'smooth binning' technique uses two SSWML calcula- 
tions with bins of two magnitudes, shifted by one magnitude 
relative to each other. This achieves the closest results to the 
input LF. 

Figure [l(] compares the LFs derives from the 3 clustered 
catalogues and the unclustered catalogue. Both the SSWML 
and 1/Knax methods derive successively steeper LFs as the 
strength of the clustering is increased, whereas the SSWML 
esimator is clearly less affected. In summary, these tests re- 
veal the superiority of the SSWML method in surveys af- 
fected by clustering. 



4.2 The Relation between the Methods 

Figure ^| reveals a interesting correspondence between the 
1/Knax and the maximum- likelihood technique. It appears 
that the smaller the SSWML bins, the closer it approximates 
the 1/Knax method. Looking back to Equation ba, the sum 
over the k objects in the survey in the numerator may be 
replaced with n pq , the number of objects in the appropri- 
ate bin. Furthermore, consider the case where there is only 
one type of galaxy. This is not a restrictive assumption, as 
one could split the surveyed galaxies by type, calculate the 
LF of each type independently prior to addition. These two 
approaches result in 



Luminosity Function Evolution by Spectral Type 13 



°r 2 



-s- 



C: 

o 



1 /V Method 

' max 

SSWML Method 



Q L J I I I I I I I I I I I I I I I I L 

-22 -20 -18 -16 -14 -12 
M — 5 log h 



1 /V 



Method 



SSWML Method 



□ 

9r 2 



3 
Q. 

o 
-e- 



Q I I I I I I I I I I I I I I I I I I I I I 

-22 -20 -18 -16 -14 -12 
M — 5 log h 



f 2 



_. 
O 

■e- 



1 /V Method 

/ max 

SSWML Method 



V 



-22 -20 -18 -16 -14 -12 
M - 5 log h 



1/V Method 

i max 

SSWML Method 



a. 2 - 



Cl - ? 



O 

■e- 



10 -18 -16 -14 -12 
M - 5 log h 



Figure 10. Luminosity functions (0.02 < z < 0.15) as a function of clustering. The top-left panel has no clustered galaxies. In the 
top-right panel, 45% of the galaxies in the B=17— 19.7 subcatalogue are in a cluster at z=0.05. In the lower panels 65% (left) and 85% 
(right) of these galaxies are clustered. 



PpQ observed 
"gal 



J2ij PijdI/dpi 



dl j dp p 



(31) 



where the sum over k has been replaced by a multiplication. 
Referring to Equation E7l the sum in the numerator may be 
replaced with n*; al , yielding 



Ppq 



prcd 

Vi 



dl/dpp 



(32) 



It appears that the iterative process is not required in this 
simplified case. The normalisation of p pq must be determined 
'/nllf, yielding 



by multiplying by n^f erved 



Ppq — 



dl j dpp 



(33) 



The final connection is Equation |22| in the limit when AM 
is small, such that fi(m) can be assumed to be constant 
across AM. As a survey has magnitude limits this can only 
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Figure 9. Clustered test catalogue of 1800 galaxies. The lumi- 
nosity function of the clustered test catalogue was determined 
using a variety of methods (see text). 



be an approximation. In this way, the inner integral may be 
approximated by a product: 



dl 

dp pq 



z p +Az/2 



v m 



-Az/2 
max AM 



n(m(z,M q ))AM^-dz 
dz 



(34) 



Substituting Equation |3j| into Equation |33| yields the famil- 
iar l/Vmax equation, 



p(z v ,M q ) = p pq 



(35) 



and we have come full circle — in the limit of small binwidth 
the SSWML and 1/Vmax methods are identical. 

Finally, we compare the results of the 1/Vmax an d SS- 
WML methods when applied to our combined redshift cat- 
alogue. Figure [ll] compares the redshift evolution of the lu- 
minosity function obtained with the 1 /V max method (equiv- 
alent to Figure 10 of Paper I, though with broader bins) and 
that determined using the SSWML method — the agreement 
is extremely good, as implied by the above discussion. In the 
rest of this paper we therefore show only the LFs estimated 
using the SSWML method, although we have checked and 
found agreement with the 1/Vmax method in every case. 
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Figure 11. Evolution of the luminosity function estimated by 
the 1/Vmax method (top) and the SSWML method (bottom). 
The luminosity function at several redshifts is depicted by the 
symbols in the legend with lcr errorbars. The long-dashed curve 
is the best-fit Schechter function determined by Loveday et al. 
(1992) without correction for Malmquist bias. All the luminosity 
function figures that follow have a similar format. 
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5 THE EVOLUTION OF THE LUMINOSITY 
FUNCTION BY SPECTRAL TYPE 

In Paper I (and Figure [H]) we demonstrated evolution in the 
shape of the overall LF since z~0.75. We present these re- 
sults again here using both the traditional 1 /Vmax technique 
and the clustering-insensitive SSWML method. A compari- 
son of these figures demonstrates that, as we argued in Pa- 
per I, the Autofib survey is not grossly affected by clustering. 
However, the SSWML techinique yields a slightly different 
and a considerably more robust estimate of the galaxy lu- 
minosity function. 

We also showed in Paper I how the steepening of the 
faint end slope with look-back time is caused primarily by 
the population of galaxies with strong [Oil] emission, which 
effectively trace the star-forming component of the field 
galaxy population. The techniques described in this paper 
allow us to directly examine the evolution of the galaxy LF 
as a function of spectral type. 



5.1 Misclassification and A-correction Error 



The analysis in §2.2 shows that the cross-correlation tech- 
nique may misclassify approximately one-fifth of the spec- 
tra in the survey by one spectral class. Therefore before 
proceeding we must check that such misclassification does 
not seriously affect our results. To address this, 20% of the 
galaxies were randomly reclassified by one class blueward 
or redward, the sampling volumes and absolute magnitudes 
re-calculated, and the luminosity functions re-determined. 
Figure |l^ presents the results of this test. These luminos- 
ity functions agree with those presented in Figure [H] to 
within the estimated errors, showing that the luminosity 
functions derived are not sensitive to the level of misclassi- 
fication present in our sample. 

To test the sensitivity to the fc-corrections themselves, 
we calculate the evolving luminosity function without con- 
verting the observed magnitudes to the rest-frames of the 
galaxies. The evolving observer's frame luminosity function 
is calculated using the SSWML method (Figure g). Here, 
the fc-corrections determined for the surveyed galaxies play 
no role in the derivation of the luminosity functions. For 
comparison, we calculate how the observed local distribu- 
tion of galaxies would appear at higher redshifts if they are 
observed in the B-band and binned by their observed B- 
band magnitude, i.e. the 'no-evolution' prediction. To per- 
form this comparison we assume a specific local mix of galax- 
ies with varying fc-corrections. As the galaxies suffer larger 
and larger fc-corrections, the bright-end cutoff luminosity 
becomes markedly fainter with redshift in both the obser- 
vations and the model. However, the model predicts fewer 
faint galaxies than are observed. This effect is analogous to 
the increase in the faint-end slope in the rest-frame luminos- 
ity functions. The number of faint galaxies must therefore 
increase with look-back time, independent of possible errors 
in the fc-corrections that we have derived. 
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Figure 12. Evolution of the luminosity function with 20% ran- 
dom re-classifications by one spectral class. 
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Figure 13. Evolution of the luminosity functions in the ob- 
server's frame. 
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5.2 Characterising the Evolution 

For convenience, we divide the sample into three broad spec- 
tral types: E/SOs, early spirals (Sab/Sbc) and late spirals 
(Scd/Sdm/starburst) in accordance with the discussion in 
§2. For each of these subsamples we compare the co-added 
spectra at low and high redshift. The spectra were coadded 
in the rest-frame with normalisations which are allowed to 
vary in order to minimise the mean dispersion between the 
various spectra. Since each spectrum covers a slightly differ- 
ent wavelength range, this is not exactly equivalent to nor- 
malising according to the total flux in each spectrum, but 
the difference only marginally affects the coadded spectra 
in the range of interest. Because of uncertainties in fluxing 
these data drawn from such diverse sources, we have made 
no attempt to flux-calibrate these data. Rather, they are in- 
tended to illustrate possible evolution of the mean spectral 
features of the galaxies. 

We also determine the LF as a function of redshift for 
each spectral type by the SSWML method, using a coarser 
sampling than for the combined sample due to the smaller 
numbers within each spectral class. As a different approach 
to characterising the LF evolution of each class, we fit a 
model for the evolution of the Schechter parameters with 
redshift using the generalised STY (SSTY) method derived 
in Equation This is done by maximising the likelihood of 
observing each subsample assuming the LF evolution takes 
the following form: 



4f(z) = <Po(i+z)^ 

L*(z) = Lt{l + Z y L 
a(z) = a + laz. 



(36) 



In the above, both the normalisation and characteristic lu- 
minosity evolve as a power of cosmic time ((f>* oc £~ 37 */ 2 
and L* oc t _37 - L/ ' 2 ), while the faint-end slope evolves linearly 
with redshift. As with a non-evolving Schechter function, the 
parameters are highly correlated. The density of galaxies of 
a given class is given by: 



cb(L,z) = 05(i + ^) 7 *~ 7i(ao+7QZ) ( il 



(37) 



exp 



L* (l + z) 



For example, if the faint-end slope does not change with 
redshift (7^=0), then <j>Q and the product aoLg jointly de- 
termine the density of faint galaxies as a function of redshift. 
The correlation of the six parameters makes the errors diffi- 
cult to interpret. However, although the trends in individual 
parameters are uncertain, those which determine the over- 
all LF shape are much better determined (at least within 
the evolutionary formalism adopted). Although the best-fit 
parameters maximise the likelihood of observing the survey, 
several sets of these parameters will produce nearly identi- 
cal trends within the absolute magnitude and redshift range 
probed. 



Furthermore, it is essential to note that the given SSTY 
model is restricted to monotonic evolution of the parameters 
with redshift. It will find the trend that the bulk of the 
galaxies of a given spectral type follow. Since the catalogue 
is numerically dominated by galaxies at moderate redshift, 
the SSTY model is most strongly sensitive to the recent 
evolution of the population. 



5.3 Early- Type Galaxies 

Galaxies with early- type (E/SO-like) spectral classifications 
represent 20% of the combined survey (331/1603 classified 
galaxies) and contribute approximately one-tenth of the vol- 
ume density of intrinsically faint galaxies and one-half of the 
galaxies with L ~ L* . 

Figure [l4| compares the mean spectra of the E/SO galax- 
ies with < z < 0.2 with those 0.2 < z < 0.5. At first sight 
the two spectra appear very similar to each other a nd also 



to the E/S spectra given in the Kennicutt atlas (Kenni- 



cutt 1992a). The most striking differences concern the 4000 



A break and [O II] emission which both appear to strengthen 
with redshift. The former effect is surprising unless there is 
a selection bias in favour of detecting galaxies with strong 
features and so cannot be considered a reliable result given 
the absence of a flux calibration. The presence of the [Oil] 
line in the coadded spectra is more interesting, although it 
could readily arise from the presence of a few later-type sys- 
tems with red spectra, as evidenced in the discussion of the 
precision of our classification system in §2.5. 

The important point here is the absence of evidence 
for recent star formation in this predominant ly spectrally- 
inactive class (e.g. a prominent H8 line, 



Zabludoff ct 



al. 1996). This indicates that the bulk of the redder sources 



in our survey have not recently undergone any star forma- 
tion on timescales of 1-2 Gyr and thus are passively evolving 
old systems over the redshift ranges probed. 

Figure ^ shows the evolution of the associated LF de- 
termined using the SSWML method in three redshift ranges, 
together with the LFs resulting from SSTY fits for the pa- 
rameters of the model described above. Neither method finds 
any convincing evolution of the E/S0 LF out to 2=0.35, 
while the small sample size at higher redshifts produces large 
uncertainties and the two methods suggest some evolution, 
but in opposite senses. The overall conclusion, from both the 
co-added spectra and the LF analysis, is that there has been 
little significant change in the properties of the early-type 
field galaxy population to at least z~0.5. 



5.4 Early- Type Spiral Galaxies 

Early-type spirals (Sab/Sbc) comprise 40% of the spectrally- 
classified galaxies (643/1603) and contribute about one-half 
of the local density of galaxies for luminosities probed by the 
combined survey. Our analysis proceeds as for the early-type 
galaxies. Figure [l6| compares the coadded spectrum for the 
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Figure 14. Coadded spectra for early-type galaxies. All the 
galaxies with E/SO spectral classifications have been coadded in 
two groups. The bold curve is the coaddition of all early-types 
with redshifts z<0.2. The light curve is for those with 0.2<z<0.5. 
Both spectra have been smoothed on a scale of 10A. 
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early spirals with z<0.2 and 0.2<z<0.5. An effect similar 
to that found in the early-types is also evident here, with 
the absorption lines in the high-redshift spectrum appear- 
ing somewhat stronger than in the low-redshift spectrum. 
The difference between the blue continuum slopes of the 
two spectra is also in the same direction as in Figure [l4| 
and presumably results from the rapid decrease in the spec- 
trograph response blueward of 3600A which preferentially 
affects the spectra of the low-redshift galaxies. Again, the 
principle result is the remarkable similarity of the diagnos- 
tic spectral features over this redshift range. We find no sig- 
nificant change in the relative contributions of the Balmer 
absorption lines or the [Oil] 3727A emission line between 
the low- and high-redshift spectra (the median EWs of [Oil] 
are similar). 



Figure 15. The SSWML luminosity functions of E/S0 galaxies 
in three redshift ranges are shown by the points. The results of 
the SSTY LF fits for z=0.1, 0.3 and 0.5 are superimposed. The 
best-fitting monotonically evolving model (found by the SSTY 
method) undcrprcdicts the number of elliptical galaxies at high 
redshift relative to the less constrained SSWML results. The Love- 
day et al. (1992) LF (dashed curve) is shown for reference. In this 
and the subsequent LF figures, we have extrapolated the SSTY 
LF fits faintward of the observed data to illustrate the trends 
more clearly. Extrapolations are defined via a dotted line. 



with the late-type spiral class discussed here and thus we 
can anticipate similar results. 

Importantly, the coadded spectra for galaxies of this 
class (Figure fl§) shows a marked change in the na ture of 



Hoj wovor, ao Figure [l7|ahowo, there i 



jdoot 



the star formation present. As originally discussed by 



lution o f the early spirals' LF with redshift , Although there 



Broad- 



hurst, Ellis & Shanks (1988), the higher redshift spectrum 



is little or no change in the number of objects with L>L* , 
the faint end of the LF steepens as we go to higher redshift. 
The absence of post-starburst features in the high redshift 
coadded spectrum indicates the bulk of the population is un- 
dergoing constant or smoothly-declining star formation on 
timescales of a few Gyr. 



5.5 Late- Type Spiral Galaxies 

As shown in Paper I, late-type galaxies appear to be the 
principal source of evolution in the shape of the overall 
galaxy LF. In Paper I, the selection method that produced 
this conclusion was based only on the EW of the [O II] emis- 
sion line. In fact this selection correlates reasonably closely 



shows much deeper HS absorption and a possible reversal of 
the Ca II line ratio suggesting the presence of a stronger He; 
these lines are signatures of more recent (~1 Gyr) star for- 
mation. Coupled with this is a strong increase in the strength 
of the [O II] emission. These diagnostics suggest such galax- 
ies are being seen at an atypical stage in their star formation 
history. 

Figure ^ reveals the [O II] evolution more clearly: at 
fixed luminosity there is an increase in the median [Oil] EW 
with redshift. This effect is superimposed upon the well- 
known tendency for fainter objects to have stronger emis- 
sion. Analysing objects in a narrow range of spectral type 
minimises the bias arising from the fact that fainter samples 
become more dominated by suc h galaxies because of the 
fc-correction (|Koo fc Kron 1992j). Taking [Oil] EW as an 
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Figure 16. Coaddcd spectra for early-type spiral galaxies. The 
bold curve is the coaddition of all early spirals with redshifts 
z<0.2, while the light curve is those with 0.2<z<0.5. Both spectra 
have been smoothed on a scale of lOA. 
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Figure 18. Coaddcd spectra for late-type spiral galaxies. The 
bold curve is the coaddition of all late spirals with redshifts 2<0.2, 
while the light curve is those with 0.2<,z<0.5. Both spectra have 
been smoothed on a scale of 10A. 
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Figure 17. SSWML luminosity function of early spiral galaxies. 
The results of the generalised STY fitting for the early-type spirals 
are superimposed for 2=0.1, 0.3 and 0.5. 
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Figure 19. Evolution of the median [Oil] EW for late-type spi- 
rals. In regions of the Mb — z where no galaxies were observed, 
a median EW of OA is plotted. 



indicator of star-formation rate, we can therefore conclude 
that in late-type spirals the star-formation rate is higher 
in lower luminosity objects at fixed redshift and higher in 
higher redshift objects at fixed luminosity. Alternatively, a 
given star-formation rate is found at higher redshift in higher 
luminosity objects. The spectra indicate the nature of star 
formation changes also, becoming more 'burst-like' in the 
distant samples. 

The evolution of the late-type spirals dominates the 
more modest changes seen in the LFs of the other types. 
Their density and faint-end slope both increase with red- 
shift, while the cutoff luminosity remains fixed. Conse- 
quently they contribute an increasing fraction at all lumi- 
nosities, progressively dominating the population at brighter 
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Figure 20. SSWML luminosity function of late spiral galaxies. 
The results of the generalised STY fitting for the late-type spirals 
are superimposed for 2=0.1, 0.3 and 0.5. 



limits and by 2~0.5 becoming the most common galaxy type 
at all luminosities below L* . 



5.6 Summary and a Check 

The evolution of the various galaxy types is summarised 
in pictorial form in Figure |2lj. This characterisation goes 
further than in the preceding sections and provides SSTY 
characterisations of the LF evolution for each of six spec- 
tral classes. With this level of subdivision the LFs are in- 
dividually highly uncertain even within the framework of 
the assumed evolutionary form (especially for the early-type 
class which is poorly represented at high redshift) . Nonethe- 
less the figure illustrates that the earlier spectral types only 
weakly evolve whereas the rate of change is greatest for the 
later types. It is also clear that the evolution takes the form 
of a steepening of the faint end of the LF with L* galaxies 
only being significantly affected at z~0.5. The parameters 
of the SSTY fits are provided in Table ^, although, in light 
of the above caveats, should clearly by used and interpreted 
with caution. 

As a final check, the SSTY LF fits for each spectral 
type are combined to give the evolution of the overall lumi- 
nosity function. This is compared to the SSWML LFs for 
the whole sample in Figure G2j and the consistency check 
demonstrates that for the large combined sample the two 
approaches yield consistent results. We conclude that the 
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Figure 21. The evolution of the luminosity function by spec- 
tral type. The curves trace schematically the best-fitting, mono- 
tonically evolving luminosity functions determined by the SSTY 
method at 2=0.1, 0.3 and 0.5. Arrows indicate increasing redshift. 



evolutionary framework assumed in the SSTY fits is suffi- 
ciently general to adequately match the observations. 



6 DISCUSSION 

The evolution of the LF by spectral type described in the 
previous section considerably extends the analyses of previ- 
ous studies, whi ch have either u sed [Oil] EW (Paper I) or 
red/blue colour (Lilly et al. 1995) to make a simple distinc- 
tion between the star-forming and quiescent galaxy popula- 
tions. Comparison of those studies with the more detailed re- 
sults presented here, however, shows broad agreement. The 
analyses based here will pave the way for more comprehen- 
sive analyses that will become possible in more extensive 
homogeneous faint spectral surveys. 

In comparing the LF evolution of the red and blue 



galaxies in their I<22 redshift survey, Lilly et al. (1995) 



found that there was very little change in the luminosity 
or number density of the red galaxies over the whole red- 
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Table 2. The evolution of the LF parameters. 



Type 
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-3 


1.56 
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-0.35 
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Sab 
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-3 


-4.44 
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0.92 


-0.99 
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Sbc 


2.80x10" 


-3 


2.89 


-19.4 
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-0.34 


Sdm / Star burst 
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-3 
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-0.97 


-1.36 


-0.79 
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we obtain for low-redshift mid- to late-type spirals in our 
SSTY parametric analysis (see Table 3), but quite consis- 
tent with the values we find for these types at z<0.5. With 
a much larger sample at lower redshifts, our results not only 
show that significant evolution is occurring at z<0.5, but 
also demonstrate the marked type-dependence, with late- 
type spirals evolving more strongly than mid-type spirals. 

Many authors have previously postulated that late-type 
galaxies are the prime movers in the evolution of the lumi- 
nosity function as appears to be manifest in the number- 
magnitude relation. Their blue colours imply that they are 
least affected by fc-corrections and thus can be obse rved to 
the largest distances for a giv en absolute magnitude 



Broad- 



hurst, Ellis & Shanks (1988) found that a luminosity func- 



tion with a constant cutoff luminosity and a faint-end slope 
which increases with redshift is consistent with both the 
slope of the observed blue numbe r co unts and the redshift 
distribution. iLacey & Silk (1991)1 and iTreyer & Silk (1994) 



both supported this finding, and conjectured that these 
changes are driven by blue, late-type galaxies which have 



disappeared since a redshift of z~0.2. Lacey fc Silk (1991) 



Figure 22. Comparison of the SSTY and SSWML luminosity 
functions. The sum of the SSTY LFs for the six classifications 
and the unclassified galaxies (the solid curves) is compared with 
the LF for all the galaxies in the survey as determined by the 
SSWML method. The Loveday et al. (1992) LF is shown as the 
short-dashed curve. Just below the solid curves are long-dashed 
curves, which show the evolution of the luminosity function with 
unclassified galaxies excluded from the analysis. The SSTY results 
are given for 2=0.1, 0.3 and 0.5. 



propose several options for their disappearance 

The excess galaxies could have produced more high- 
mass stars proportionally than modern galaxies. This bias 
toward high-mass stars would make them brighter while 
stars were forming, but the galaxies would quickly fade after 
the high-mass stars became supernova, leaving few stars to 
be observed today. An alternative explanation is that these 
galaxies have mer ged into the galaxie s that we observe today. 
White (1990) and Efstathiou (1990) argue that such an evo- 
lutionary history is a natural consequence of the turnaround 
of larger and larger mass scale with time in a hierarchical 
universe. However, large amounts of merging wou ld leave 
traces in toda y's galaxies which are not observed (Toth & 



shift ra |nge 0<z<l, a point which we can confirm both in Ostriker 1992j ). Furthermore, these excess galaxies are un 



terms of the derived LF and the mean spectral cha racter- 
istics. Amongst the blue galaxies Lilly et al. (1995) found 
substantial evolution at z>0.5, so that by 0.5<z<0.75 the 
blue LF had brightened by about 1 mag. They were non- 
committal about the evolution of the blue LF at lower red- 
shifts, due to the small number of objects at z<0.2 (35) 
with which to compare their 0.2<z<0.5 LF. Nonetheless, 
their Schechter function parameters indicated that the blue 
galaxies with a<0.5 (which are mostly at the high end of our 
redshift range) have a LF with a faint-end slope a between 
— 1.3 and —1.6. This is as steep or steeper than the value 



likely to be the progenitors of today's galaxies as they are 
more weakly clustere d and more dense than galaxies today 
(Babul & Rees 1991). However, there is no convincing evi- 



dence in either our early-type spectra or L Fs for a significan t 
build up of more massive early types (c.f. Dalcanton(1993)) 
as would be the case in the merger hypotheses. 

Babul & Rees (1991) propo se a further explanation (elu- 
cidated also in Efstathiou 1992). Although small mass haloes 
(~ 1O 9 M0) collapse and virialise before the nascent haloes of 
L* galaxies, the UV flux produced by quasars may keep the 
gas in small haloes ionised until Only then can these 
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small galaxies begin to form stars. Such stars would form 
quickly as in a starburst galaxy. During the starburst these 
galaxies would appear irregular and their spectra would be 
similar to those of present-day late-type spirals (i.e. domi- 
nated by short-term stellar activity) . Supernovae would blow 
out the gas, halting the star formation rate, and the galaxy 
would begin t o fade (more quickly in the B-band than in 
the JC-band). Babul & Rees (1991) also proposed that in 



low-pressure environments the gas may escape the galaxy 
entirely, while in intermediate and high-pressure regions, 
some gas might return to the galaxy, possibly fuelling further 
events. The few galaxies that could still be observed today 
would be in the high-pressure regions clustered near more 
luminous galaxies. Gigayears of phase mixing could trans- 
form these irregular starburst galaxies into today's popula- 
tion of dwarf elliptical galaxies. However, the vast majority 
of these remnant irregulars would fade below detection lim- 
its. Although the clear proof of this picture would be the 
identification of an abundant population of local remnants, 
presumably via deep infrared imaging, the short-term na- 
ture of star formation in the coadded spectra of the domi- 
nant evolutionary component provides some support for this 
picture. With the larger databases anticipated in the forth- 
coming deep surveys, it will hopefully become possible to 
tie the timescales observed in the redshift evolution to those 
indicated by the detailed spectral features in the coadded 
spectra. 



7 CONCLUSIONS 

This paper investigates the evolution of the galaxy luminos- 
ity function for galaxies of different spectral types. The anal- 
ysis is based on the Autofib redshift survey described in 



Ellis 



et al. (1996) and uses new techniques developed for spectral 



classification and for recovering the luminosity function in 
the presence of both evolution and clustering. 

We classify the redshift survey galaxies into six spec- 
tral types using a cross-correlation method that simulations 
show correctly assigns the type for 80% of the sample and is 
out by at most one type more than 90% of the time. These 
types compare well with those derived from the galaxies' 
colours and, where available, with Hubble Space Telescope 
morphologies. 

We describe an extension of the step-wise maximum 
likelihood method for estimating luminosity functions, rig- 
orously testing the method and comparing its performance 
to the standard 1/Vmax method. Significant advantages to 
the new SSWML method are found in reducing the sensitiv- 
ity of t he results to clustering. We also desc ribe an extension 



to the 



Sandage, Tammann & Yahil (1979) 



method that en- 



ables us to recover the best fit to a parametric model for the 
luminosity function including its evolutionary behaviour. 

Applying these methods to the Autofib redshift survey 
we obtain the following main results: 

1. There is no significant evolution in the mean spec- 



trum or luminosity function of early-type galaxies out to 
at least z~0.5. This significantly constrains the continued 
production of such galaxies (e.g. via mergers of star-forming 
galaxies over the last few Gyr). 

2. Early-type spiral galaxies show relatively modest evo- 
lution which is well-characterised by a simple steepening of 
the faint end of their luminosity function with redshift. Early 
spirals show little evidence for evolution of either L* or (j>* 
over this redshift range. Their spectra at all redshifts sam- 
pled are consistent with smooth changes in the overall star 
formation history. 

3. Out to z~0.5, the overall evolution of the galaxy 
population is dominated by the evolution of late- type spirals. 
Their luminosity function not only steepens and brightens 
(evolution of both a and L*) but there are also signs of 
significant density evolution (a rapid increase in (/>*). There 
appears to be a qualitative change in the spectra of late 
type spirals with redshift. In addition to a rise in the median 
[Oil] equivalent width, a greater proportion of high redshift 
sources show Balmer line ratios indicative of short-term star 
formation. Coupled to the increase in star-formation rate 
at lower luminosities at fixed redshift, this implies that a 
given star-formation rate is found at higher redshift in higher 
luminosity objects. 
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